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We review the energy spectrum and transport properties of several types of one- 
dimensional superlattices (SLs) on single-layer and bilayer graphene. In single-layer 
graphene, for certain SL parameters an electron beam incident on a SL is highly 
coUimated. On the other hand there are extra Dirac points generated for other SL 
parameters. Using rectangular barriers allows us to find analytic expressions for 
the location of new Dirac points in the spectrum and for the renormalization of the 
electron velocities. The influence of these extra Dirac points on the conductivity 
is investigated. In the limit of 5-function barriers, the transmission T through, 
conductance G of a finite number of barriers as well as the energy spectra of SLs 
are periodic functions of the dimensionless strength P of the barriers, P5{x) — 
V{x)/hvp, with V F the Fermi velocity. For a Kronig- Penney SL with alternating sign 
of the height of the barriers the Dirac point becomes a Dirac line for P = 7r/2 -I- mr 
with n an integer. In bilayer graphene, with an appropriate bias applied to the 
barriers and wells, we show that several new types of SLs are produced and two of 
them are similar to type I and type II semiconductor SLs. Similar as in single-layer 
graphene extra "Dirac" points are found. Non-ballistic transport is also considered. 
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1. Introduction 



Since the experimental realisation of graphene ( Novoselov et al. 2004 ) in 2004 



this one-atom thick layer of carbon atoms has attracted the attention of the scien- 
tific world. This attraction pole was created by the prediction that the carriers in 
graphene behave as massless relativistic fermions moving in two dimensions. The 
latter particles, which are described by the Dirac- Weyl Hamiltonian, possess inter- 
esting properties such as a gapless and linear-in-wave vector electronic spectrum, 
a perfect transmission, at normal incidence, through any potential barrier, i.e., the 



Klein paradox (Klein 1929 Katsnelson et al. 2006 Pereira Jr et al. 2010 Roslyak 



et al. 2010), which was recently addressed experimentally (Young & Kim 2009 



2007 



Huard et al. 2007), the zitterbewegung (Schliemann et al 



Zawadzki 2005), etc., see Ref. (Castro Neto et al. 2009) and (Abergel et al. 



2005 Winkler et al. 



2010 ) for recent reviews. On the other hand, in bilayer graphene the carriers exhibit 
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a very different but extraordinary electronic betiaviour, suclr as being cliiral (Mc- 



Cann 2006 Katsnelson et al. 2006) but with a different pseudospin (=1) than in 



single-layer graphene (=1/2). Although their spectrum is parabolic in wave vector 
and also gapless, it is possible to create an energy gap by applying a perpendicular 
electric field on a bilayer graphene sample ( Castro et al. 2007 1 . This allows one to 



200761 



electrostatically create quantum dots in bilayer graphene ( Pereira Jr et al. 
and enrich its technological capabilities. 

In previous work we studied the band structure and other properties of single- 



layer and bilayer graphene (Barbier et al. 2008 2009 &) in the presence of one- 
dimensional (ID) periodic potential, i.e., a superlattice (SL). SLs are known to 
be useful in altering the band structure of materials and thereby broadening their 
technological applicability. 

The already peculiar, cone-shaped band structure of single-layer graphene can be 
drastically changed in a SL. An interesting feature is that for certain SL parameters 



the carriers are restricted to move along one direction, i.e. they are coUimated (Park 



et al. . 2009a ) . Furthermore, it was found that for other parameters of the SL instead 



of the single-valley ( the K or iiT'-point) Dirac cone, "extra Dirac points" appeared 
at the Fermi level in addition to the original one ( Ho et al. 2009 ) . The latter extra 



Dirac points are interesting because of their accompagning zero modes (Brey & 



Fertig 20091 and their influence on many physical properties such as the density 



of states (Ho et al. 20091, the conductivity (Barbier et al. 2010 Wang & Zhu 



2010), and the Landau levels upon applying a magnetic field (Park et al. 2009 & 



Sun et al. 2010). 



One can also obtain "extra Dirac points" in bilayer graphene SLs. The possibility 



of locally altering the gap ( Castro et al. 2007 ) of bilayer graphene by applying a 



bias is another way of tuning the band structure. In this review we classify these 
SLs in four types. Another interesting result of applying a bias locally is that sign 



flips of the bias introduce bound states along the interfaces (Martin et al. 2008 



Martinez et al. 2009 ) . These bound states break the time reversal symmetry and 



are distinct for the two K and K' valleys; this opens up perspectives for valley-filter 



devices ( San- Jose et al. 2009 1 



In this review we will use the following methods to describe our findings. For 
both single-layer and bilayer graphene we will use the nearest neighbour, tight- 
binding Hamiltonian in the continuum approximation, and restrict ourselves to 
the electronic structure in the neighbourhood of the K point. We then apply the 
transfer-matrix method to study the spectrum of and transmission through various 
potential barrier structures, which we approximate by piecewise constant potentials. 
We consider structures with a finite number of barriers and SLs. 

We will study ballistic transport in systems with a finite number of barriers using 
the two-probe Landauer conductance while in a SL (infinite number of barriers) we 
will evaluate the spectrum and the diffusive conductivity, i.e., we will study non- 
ballistic transport. 

The work is organized as follows. In Sec. [2] we investigate various aspects of 
ballistic transport through a finite number of barriers on single-layer graphene as 
well as the spectrum of SLs, with emphasis on collimation and extra Dirac points 
and their infiuence on non-ballistic transport. In Sec. [3] we carry on the same stud- 
ies, whenever possible, on bilayer graphene. In addition, we consider various types 
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of band alignments in the presence of a bias that can lead to different types of 
heterostructures and SLs. We make a summary and concluding remarks in Sec. |4j 

2. Single-layer graphene 

We describe the electronic structure of an infinitely large, flat graphene flake by 
the nearest-neighbour tight-binding model and consider wave vectors close to the 
K point. The relevant Hamiltonian in the continuum approximation is "H — vpff ■ 
■p + V^ + mv^dzj with p the momentum operator, V the potential, 1 the 2x2 unit 
matrix, a — {ax<Jy), (j^ the Pauli-matrices and vp ~ lO^m/s the Fermi velocity. 
Explicitly "H is given by 



U 



V + mvp —ivph{dx-'idy) 
—ivFh{dx + idy) V — mvp 



(2.1) 



The mass term is in principle zero in the nearest-neighbour, tight-binding model 



but due to interaction with a substrate ( Giovannetti et al. 2007 1 an effective mass 



term can be induced and results in the opening of an energy gap. Recently there 
have been proposals to induce an energy gap in single-layer graphene, and it is 
appropriate that we consider this mass term where relevant. In the presence of 
a ID rectangular potential V{x), such as the one shown in Fig. [l| the equation 
{H — E)^p = admits (right- and left-travelling) plane wave solutions of the form 
?/'i,r(a:)e*'^f^ with 



Ai^) - [x + iky) '^'(^^ 



-A + iky 



(2.2) 



here A = [(e — u{x))'^ 



i^]^/^ is the X component of the wave vector, e ~ 



EL/hvp, uix) — V{x)L/hvp, and fi — mvpL/h. The dimensionless parameters e, 
u{x) and fj, scale with the characteristic length L of the potential barrier structure. 
For the single or double barrier system this L will be equal to the barrier width 



while for a SL it will be its period. Neglecting the mass term one rewrites Eq. (2.2) 
in the simpler form 



1 



V'i(a;) 



1 

-se' 



with A = [(£ — u{x))^ — kyY^^, tan0 = ky/\, and s — sgn(e — u{x)). 



(2.3) 



(a) A single or double barrier 

The model barriers and wells we consider are shown in Fig.[l|a). It is interesting 
to look at the tunneling through such barriers, which was previously studied by 



(Katsnelson et al. 20061 for a single barrier. This was later extended to massive 



2008) 



electrons with spatially varying mass (Gomes & Peres 

Transmission. To find the transmission T through a square-barrier structure 
one first observes that the wave function in the jth region ipjix) of the constant 
potential Vj is given by a superposition of the eigenstates given by Eq. (2.2), 



ipjix) = Aj^pr■j + Bji>u 



(2.4) 
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A V(x) 



"b , 



(a) 



(b) 



Figure 1. (a) A ID potential barrier of height Vb and width Wb- (b) A single unit of a 
potential well next to a potential barrier. 

The wave function should be continuous at the interfaces. This boundary condition 
gives the transfer matrix Mj relating the coefficients Aj and Bj of region j with 
those of the region j + 1 in the manner 



(2.5) 



By employing the transfer matrix at each potential step we obtain, after n steps, 
the relation 

TT ^r (An 



In the region to the left of the barrier we assume Aq — \ and denote by = r the 
reflection amplitude. Likewise, to the right of the nth barrier we have i?„ = and 
denote by = t the transmission amplitude. 

The transmission probability T can be expressed as the ratio of the transmitted 
current density jx over the incident one, where jx = This results in 

T — (A'/A)|ip, with A'/A the ratio between the wave vector A' to the right and A 
to the left of the barrier. If the potential to the right and left of the barrier is the 
same we have A' — A. For a single barrier the transmission amplitude is given by 
T = |ip = lA^iil"^, with Nij the elements of the transfer matrix J\f. Explicitly, t 
can be written as 

l/t = cos(A6M^6) - iQ siii{XbWb), 
Q = {eoEb - ky - /io/i6)/AoAb; 

the indices and h refer, respectively, to the region outside and inside the barrier 
and £{, = e — M. A contour plot of the transmission is shown in Fig. [2ja). We 
clearly see: 1) T = 1 for = which is the well-known Klein tunneling, and 2) 
strong resonances, in particular for _B < 0, when XhWh — nn, which describe hole 
scattering above a potential well. 

In the limit of a very thin and high barrier, one can model it by a (5-function 



barrier V{x)/hvF = P5{x). Using Eq. (2.7| for t gives (Barbier et at. 2009a) 

T = + sin2ptan2 0], (2.8) 

with tan0 — ky/X^ the angle of incidence. Notice that this transmission is inde- 
pendent of the energy and is a periodic function of P. The latter is very different 
from the non-relativistic case where T is a decreasing function of P. A contour plot 
of the transmission is shown in Fig. [2|b) and T = \ for (/) « which is nothing else 
than Klein tunneling. Notice also the symmetry r(7r — P) = T{P). 
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Figure 2. (a) Contour plot of the transmission through a single barrier with n = 0, Wb = L, 
and itt, — 10. (b) As in (a) for a single (5-function barrier with fi = and u{x) = PS{x); 
the transmission is independent of the energy, (c) As in (a) for two barriers with fj, — 0, 
Ub = 10, it„ = 0, Wb = 0.5L, and = L. (d) Spectrum of the bound states vs ky for 
a single (L = 1, solid red line), two parallel (dashed blue curves), and two anti-parallel 
(green dash-dotted curves) (5-function barriers (L is the inter- barrier distance). 

For two barriers the system becomes a resonant structure, for which it was found 
that the resonances in the transmission depend mostly on the width Ww of the well 



between the barriers (Pereira Jr et al. 2007a I. A plot of the transmission is shown 
in Fig. [2][c) . In the limit of two parallel (5-function barriers of equal strength P we 
obtain the transmission 



T = [l + tan^ 0(cos Aq sin 2P ~ 2s sin Ao sin^ P/ cos ( 



(2.9) 



The case of two anti-parallel (5-function barriers of equal strength is also interesting. 
The relevant transmission is 



T = [ cos^ Ao + sin^ Aq ( 1 - sin^ (^ cos 2P) V cos'' . 
Conductance. The two-terminal conductance is given by 



Go / T{Ef,' 

I-tt/2 



(2.10) 



(2.11) 



with Go = iEpLyC^ / {vph"^) for single-layer graphene, and Ly the width of the 
system. For a single and double barrier, the transmission through which is plotted 
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in Fig.[2ja) and[2]jc), the conductance G is shown in Fig.jsjjb) and exhibits multiple 
resonances despite the integration over the angle <f>. 

Taking the limit of a (5-function barrier leads to G periodic in P and given by 



G/Go = 2 [l - artanh(cos P) sin P tan P] / cos^ P. 
For one period G is shown in Fig. [sf^a) . 



(2.12) 





10 

EL / hv^ 



Figure 3. (a) Conductance G vs strength P of a 5-function barrier in single-layer graphene; 
the conductance is independent of the energy, (b) Conductance G vs energy for the single 
(solid blue curve) and double (dashed green curve) square barrier of Fig. [2]^a) and[2|c). 



Bound states. For ki 



Mo 



> the wave function outside the barrier (well) 



becomes an exponentially decaying function of x, i){x) oc exp{±|/ejx]_with_[fc^ 



[fc^-|- /iQ — g^]-'-/^. Localized states form near the barrier boundaries (Pereira Jr et al. 



2006); however, they are propagating freely along the y-direction. i'he spectrum o: 



these bound states can be found by setting the determinant of the transfer matrix 
equal to zero. For a single potential barrier (well) it is given by the solution of the 
transcendental equation 



|Ao|Afc cos(A{,VFb) + {kl + ^oA^fc - £(e - ")) sin(AbWh) = 0. 



(2.13) 



In Fig. |4](b) these bound states are shown, as a function of ky^ by the dashed blue 
(red) curves. 

An interesting structure to study is that of a potential barrier next to a well but 
with average potential equal to zero, considered by (Arovas et al. 2010). This is 



the unit cell (shown in Fig. [T|b)) of the SL we will use in Sec. [3| where extra Dirac 
points will be found. In FigTWa) the Dirac cone outside the barrier is shown as a 
grey area, inside this region there are no bound states. Superimposed are grey lines 
corresponding to the edges of the Dirac cones inside the well and barrier that divide 
the (i?, ky) plane into four regions. Region I corresponds to propagating states inside 
both the barrier and well while region II (III) corresponds to propagating states only 
inside the well (barrier). In region IV no propagating modes are possible, neither 
in the barrier nor in the well. For high thin barriers, region I will become a thin 
area adjacent to the upper cone, converging to the dark green line in the limit of 
a i5-function barrier. Figure |4][b) shows that the bound states of this structure are 
composed of the ones of a single barrier and those of a single well. Anticrossings 
take place where the bands otherwise would cross. The resulting spectrum is clearly 
a starter of the spectrum of a SL shown in Fig. |4|^d). 
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Figure 4. (a) Four different regions for a single unit of Fig. [l]^b) with ui, — 24, u^, — 16, 
Wb = 0.4 and W-w = 0.6. The green line corresponds to region I in the limit of a 5-function 
barrier, (b) Bound states for a single barrier (dashed blue curves) and well (dashed red 
curves) and the combined barrier-well unit (black curves), (c) Contour plot of the trans- 
mission through a unit with ^ — 2, Ub = —u^ = 20 and Wb ~ Ww = 0.5; the red curves 
show the bound states, (d) Spectrum of a SL whose unit cell is shown in Fig. 1(b), for 
kx = Q (blue curves) and k^L = 7r/2 (red curves). 



In the limit of (5-function barriers and wells the expressions for the dispersion 
relation are strongly simplified by setting /i = in all regions. For a single J-function 
barrier the bound state is given by 

£ = sgn(sinF)|fcy|cosP, (2.14) 

which is a straight line with a reduced group velocity Vy] the result is shown in 
Fig. ^d) by the red curve. Comparing with the single-barrier case we notice that 
due to the periodicity in P, the ^-function barrier can act as a barrier or as a well 
depending on the value of P. 

For two (5-function barriers there are two important cases: the parallel and the 
anti-parallel case. For parallel barriers one finds an implicit equation for the energy 

|A'cosP + esinP| = |e~^'/cy sinP|, (2.15) 

where A' — |Ao|, while for anti-parallel barriers one obtains 

fc2sin2p = A'V(l-e-2^'). (2.16) 
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For two (aiiti-)parallel ^-function barriers we have, for each fixed ky and P, two 
energy values ±e, and therefore two bound states. In both cases, for P = wk the 
spectrum is simplified to the one in the absence of any potential e = ±|A:j,|. In 
Fig. |2jd) the bound states for double (anti-)parallel J-function barriers are shown, 
as a function of kyL, by the dashed blue (dash-dotted green) curves. For anti- 
parallel barriers we see that there is a symmetry around = 0, which is absent 
when the barriers are parallel. 



(6) Superlattice 

Now we will consider the system of a superlattice with a corresponding ID 
periodic potential, with square barriers, given by 

oo 

V{x) ^Vo J2 - •?'^) - ©(^ - -i^ - ^b)]- (2.17) 

j=-oo 

with Q{x) the step function. The corresponding wave function is a Bloch function 
and satisfies the periodicity condition ipiL) = '0(0) ^^^{ikx), with k^ now the Bloch 
phase. Using this relation together with the transfer matrix for a single unit i^{L) = 
A4ip{0) leads to the condition 

det[7W - exp(ifc^)] = 0. (2.18) 
This gives the transcendental equation 

cosfc^r — cos XwWw cos XbWb — Q sin XwWw sin XbWb, (2-19) 



from which we obtain the energy spectrum of the system. In Eq. ( 2.19 1 we used the 
following notation: 

eu]=e + uWb, eb = e-uWw, u = VoL/hvF, Wb^w -^Wb,w/L, 
Au, = [e^ - - M™]^^^) Ab = [e^ - fc^ - ^ft]^/^, Q ^ {sweb - ky - ^bfJ-w)/Xu]Xb. 

Numerical results for the dispersion relation E{ky) are shown in Fig. |4jd). We 
see the appearance of bands (green areas) which for large ky values collapse into 
the bound states (where the red and blue curves meet) while the charge carriers 
move freely along the y direction. 



(c) Collimation and extra Dirac points 



As shown by various studies, carriers in graphene SLs exhibit several interesting 
pecularities that result from the particular electronic SL band structure. In a ID SL 



it was found that the spectrum can be altered anisotropically (Park et al.. 2008a 



Bliokh et al. 2009). Moreover, this anisotropy can be made very large such that 



for a broad region in k space the spectrum is dispersionless in one direction, and 
thus electrons are collimated along the other direction (Park et al. 2009a). Even 



more intriguing was the ability to split off "extra Dirac points" (Ho et al. 20091 



with accompanying zero modes ( Brey & Fertig 2009 ) which move away from the 
K point along the ky direction with increasing potential strength. Here we will 
describe these phenomena for a SL of square potential barriers. 
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We start by describing the collimation as done by (Park et al. 2009a); sub- 
sequently we will find the conditions on the parameters of the SL for which a 
collimation appears. It turns out that they are the same as those needed to create 
a pair of extra Dirac points. 




\\ 
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Figure 5. The lowest conduction band of the spectrum of graphene near the K point in 
the absence of SL potential (a), (b) and in its presence (c), (d) with u — in. (a) and (c) 
are contour plots of the conduction band with a contour step of 0.5 hvp/L. (b) and (d) 
show slices along constant kyL = 0, 0.2, QA-k. 



Following ( Park et al. 



2009 a 



the condition for collimation to occur is /g^ ^ 



0, where the function a{x) = 2 u{x')dx' embodies the influence of the potential, 
s = sign(£) and s — s\gii{kx). For a symmetric rectangular lattice this corresponds 
to it/4 = UTT. The spectrum for the lowest energy bands is then given by (Park 



et al. 20086) 



with // being the coefficients of the Fourier expansion e'"*^^^ = z^;=-oo ■ 
The coefficients /; depend on the potential profile V{x), with |/;| < 1. For a sym- 
metric SL of square barriers we have // = usin(/7r/2 — u/2)/{Pv? — The 
inequality |/;| < 1 implies a group velocity in the y direction Vy < vp which can be 



(2.20) 

Jl^ 



seen from Eq. (2.201. 

In Fig. [5jb),(d) we show the dispersion relation E vs for u = 47r at constant 
ky. As can be seen, when a SL is present in most of the Brillouin zone the spectrum, 
partially shown in (c), is nearly independent of ky. That is, we have collimation 
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of an electron beam along the SL axis. The condition u — V^L/hvp = Am: shows 
that altering the period of the SL or the potential height of the barriers is sufficient 
to produce collimation. This makes a SL a versatile tool for tuning the spectrum. 
Comparing with Figs.[5|a), (b) we see that the cone-shaped spectrum for m = 0, is 



transformed into a wedge-shaped spectrum (Park et al.\ 2009a 



We will compare this result now with an other approximate result for the spec- 
trum, where we suppose e small instead of ky small. We start with the transcenden- 
tal equation (2.19). As we are interested in an analytical approximate expression 



for the spectrum, we choose to expand the dispersion relation around e = up to 
second order in e. The resulting spectrum is 



e±=± 



4|a2|2 [fc2sin2(a/2) 



2 • 2 

a sm 



(fc./2)] 



/c^asina + o^u-^/lQ - 2klu'^ sin^(a/2) 



1/2 



(2.21) 



with a = [v? /A — ky]^/^. In order to compare this spectrum with that by (Park 



et al. 2009a), we expand Eq. (2.191 for small k and e; this leads to 



(2.22) 



This spectrum has the form of an anisotropic cone and corresponds to that of 
Eq. (2.20) for I = (higher I correspond to higher energy bands). In Fig. [6](a), 
(b) we see that the cone-shaped spectrum in (a), for u = 0, is transformed into 
a anisotropic spectrum in (b), for u = 4.57r, having peculiar extra Dirac points. 
These extra Dirac points cannot be described by a spectrum having an anisotropic 
cone-shape, therefore we compare the two approximate spectra. In Fig.[6](c), (d) we 
show how Eq. (2.21)) and Eq. (2.22) differ from the "exact" numerically obtained 



spectrum. From this figure one can see that Eq. (2.21) describes the lowest bands 
rather well for e < 1, while Eq. (2.22) is sufficient to describe the spectrum near 



the Dirac point. The former equation will be useful! when describing the spectrum 
near the extra Dirac points and we will use it to obtain the velocity. 

We now move on to another important feature of the spectrum, the extra Dirac 
points first obtained by (Ho et al. 2009) using tight-binding calculations. These 



extra Dirac points are found as the zero-energy solutions of the dispersion relation 
in Eq. (2.19) for zero energy (Barbier et al. 20101. 

In order to find the location of the Dirac points we assume fc^, = 0, e = 0, 
fJ'b — fJ-w — 0, and consider the special case of Wt = = 1/2 in Eq. (2.19). The 
resulting equation 



1 



■ A/2 + [(«V4 + fc^)/(«V4 - kl)] sin' A/2, 



(2.23) 



has solutions for v? / A — ky — v? / A^ 
oi ky = (at the Dirac points) as 



- ky or sin^ A/2 = 0. This determines the values 



= ± 



4j2^2. 



(2.24) 



the extra Dirac points are for j ^ 0. For a SL spectrum symmetric around zero 
energy, the extra Dirac points are at £ = 0. We expect from the considerations of 
Sec.[2]|b]) (and Fig.|4](b)) that for unequal barrier and well widths this will no longer 
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be true. Indeed, in such a case the extra Dirac points shift in energy, as seen in 
Figjlj^d), and their position in the spectrum is given, for kx = 0, by (Barbier et al. 
20l0f 

2 



e,,„ - -(1 - 2W,) + ^[^^- ^2 



± 



1/2 



(2.25) 



where j and m are integers, and m ^ corresponds to higher and lower crossing 



points. Also, perturbing the potential with an asymmetric term, as done by (Park 

(b) 



et al.. 20096), leads to qualitatively similar results. 

(a) 




Figure 6. The spectrum of graphene near the K point in the absence of a SL (a) and in 
its presence (b) with u — 4.57r. (c) and (d) The SL spectrum with u = lOvr, the lowest 
conduction bands are coloured in cyan, red, and green for, respectively, the exact, and the 



approximations given by (c) Eq. (2.211 and (d) Eq. (2.221, respectively. The approximate 
spectra are delimited by the dashed curves. 



An investigation of the group velocity near the (extra) Dirac points is appropri- 
ate for understanding the transport of carriers in the energy bands close to zero en- 
ergy. Near the extra Dirac points the group velocity tends to renormalise differently 
as compared to the original Dirac point. Near them v is oriented along the y direc- 



tion, while near the latter one v is oriented along the x direction (Ho et al.\ 2009 j ) 



The group velocity near the extra Dirac points can be calculated from Eq. (2.21 ). At 
the jth extra Dirac point the magnitude of the velocity v/vp = {de/dk^, de/dky) 
is given by 

Vx/vF = IGn^f cos{kx/2)/u^ 
Vy/vF = {uy4-^4f7T^)/u', 
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while at the main Dirac point it is given by Vx/vp = 1 and Vy/vp = 4sin(M/4)/u. 
The dependence of the velocity components on the strength of the potential barriers 
is shown in Fig. [7| From this figure we observe that new extra Dirac points emerge 
upon increasing u = VqL/Hvp (consistent with Eq. ( 2.24| )) and decreases while 



Vy increases. The Dirac point itself, however, shows a different behaviour upon 
increasing u, namely Vx = vp constant, and Vy is here a globally decaying function 
showing Vy = for periodic values of u, u ^ Anir, with n a nonzero positive integer. 



0.5 




Figure 7. The group velocity components Vy and at the Dirac point j = (shown, 
respectively, by the solid blue and the dot-dot-dashed red curve), and at the extra Dirac 
points j = 1, 2, 3 (shown, respectively, by the dot-dashed blue and the dashed red curves) 
as a function of the barrier parameter u = VoL/hvp. 



Conductivity. We now turn to the transport properties of a SL and look at 
the influence of these extra Dirac points on the conductivity. The diffusive dc con- 
ductivity for the SL system can be readily calculated from the spectrum if we 



assume a nearly constant relaxation time t{Ep) = Tp. It is given by (Charbonneau 



et al. 1982) 



<Jpiu{Ep) = 'LJ^IL ^ u„^u„^/„k(l - /nk), (2.27) 

n,k 

with A the area of the system, n the energy band index, ^, v = x, y, and /„k = 
l/[exp{f3{Ep — Enk)) + 1] the equilibrium Fermi-Dirac distribution function; /3 = 
l/fcsT and the temperature enters the results through the dimensionless value for 
/3 which is /3 = hvp /ksTL = 20. 

For comparison we first look at the conductivity tensor at zero temperature and 
in the absence of a SL. For single-layer graphene the conductivity is given by 

<J^,^,iep)/ao=ep/ATr (2.28) 

with (Jo e^/^j 

In Figs.jsf^a), (b) the conductivities axx and ayy are shown for a SL as functions 
of the energy. Notice that for small energies the slope of the conductivity ayy is 
tunable to a large extent by altering the parameter u of the SL. The dashed blue 
curves correspond to u = 47r and the rather flat dispersion in the y direction for 
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the lowest conduction band (see Fig. [5]jc,d)) translates to a small ayy (for energies 
EL/hvp < 1) compared to the conductivity in the absence of a SL. The solid 
red curves on the other hand correspond to u = 67r and due to the extra Dirac 



points, which have a rather flat dispersion in the x direction (Ho et al. 2009), the 
conductivity ayy is large. 
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Figure 8. (Color online) Conductivities, a^^ in (a) and ayy in (b), vs Fermi energy for a 
SL on single-layer graphene with u = An and 6tt for, respectively, the blue dashed and 
red solid curves. In both cases Wb = Ww = 0.5. The dash-dotted black curves show the 
conductivities in the absence of the SL potential, a^^ = ayy = epo'o/^'K. 



(d) Dirac lines 




Figure 9. (a) Schematics of Kronig-Penney SL on single-layer graphene. (b) Extended 

Kronig-Penney SL. 

In an effort to simplify the expressions for the dispersion relation we replace, as 
we did for the few-barrier structures, the SL barriers by (5-function barriers. The 
square SL potential is then approximated by 

oo 

V{x) = P ^ S{x-jL). (2.29) 

j = -oo 

This potential leads to the dispersion relation 

coskx = cos A cos P + (e/ A) sin A sin P, (2.30) 

which is periodic in P. This is in sharp contrast with that for standard electrons 
which is not periodic in P and which in our notation reads 

cos k-^ = cos A' + ii-iP/X') sin A', (2.31) 
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where fi = mvpL/h and A' = [2^£ — ky]-^/^. As can be seen from Fig. 10 a), the 
energy band near the Dirac point has the interesting property that it becomes nearly 
flat in kx, forming a plane, for large ky. The angle which the asymptotic plane makes 
with the zero-energy plane depends on P and the group velocity Vy corresponding 
to this asymptotic plane varies from —vp to vp in each period m: < P < (n + 1)it. 
Notice that no extra Dirac points are found and the reason is the same as that 
for the asymmetric SL potential, i.e., the extra Dirac points shift away from zero 
energy. Alternatively, we can try to shed some light by comparing with Sec. [2}|b]), 
where it is explained that the bound states for a single unit of the SL potential are 
similar to those of the combined single barrier and well. In the region where the 
bound states cross (denoted by I in Fig.|4][a)) anti-crossings occur and corresponding 
crossings in the SL spectrum (extra Dirac points) are expected. In the limit of a 
(5-function barrier this region is reduced to a line (the dark green line in Fig. |4][a)). 
This prevents anti-crossings from occurring and in this way no extra Dirac points 
are expected. 




Figure 10. (a) Spectrum for a Kronig-Penney SL with P = OAir. The blue and red curves 
show, respectively, the fcj; = and kx — it/L results which delimit the energy bands (green 
coloured regions), (b) Spectrum for an extended Kronig-Penney SL with P = n/2. Notice 
that the Dirac point has become a Dirac line. 



Extended Kronig-Penney model. To re-establish the symmetry between electrons 
and holes, as in the case of square barriers with Wb = Ww, we can use alternating- 
in-sign (5-function barriers. The unit cell of the periodic potential contains one such 
barrier up, at a; = 0, followed by a barrier down, at x = L/2, see Fig. [9f^b). The 
potential is given by 

C30 

j=-oo 

and is the asymptotic limit of the potential shown in Fig.[l|b). The resulting transfer 
matrix leads to the dispersion relation 

cos kx = cos A - (2fc^/A2) sin2(A/2) sin^ P. (2.33) 

This dispersion relation is periodic in P. As shown in Fig. |10[ b) no extra Dirac 
points occur, but for the particular case of P = (n + l/2)7r, n an integer, the 



Article submitted to Royal Society 



Review: Single-layer and bilayer graphene SLs 



15 



spectrum shows an interesting feature: for all ky we see that Eq. ( 2.33 1 has a solution 



with e — kx — 0, which means the Dirac point at = ky = turned into a Dirac 
line along the ky axis. If we take ky not too large (of the order of fc^;), this spectrum 
has a wedge structure as was also found for rectangular SLs. For ky oo, though, 
the spectrum becomes a horizontal plane situated at e = 0. We can generalize this 
model by taking the distance W between the two barriers of the unit cell not equal 
to L/2. This was done by (Ramizani Masir et ai, 2010, unpublished work). They 
found an approximate analytic expression for the dispersion given by 

e K [kl + Fkl]^/"^, F ^W^ + {L~Wf + 2WiL~W)cos{2P). (2.34) 

This dispersion has the shape of an anisotropic cone with a renormalized velocity 
in the y direction. Comparing with Eqs. (2.20) and (2.22), we observe that the 
condition for coUimation and the velocity renormalization in the y direction is 
quite different for square barriers. For instance, in the extended KP model, with 
W — L/2, we find Vy/vp = | cosP| while for square barriers the result is Vy/vp ~ 
sm{u / 4:) / {u / 4) . The latter means that if we consider P = u/4, the velocity in the y 
direction is maximum Vy = vp for P = (l/2 + 7i)7r in the extended KP model while 
for square barriers Vy = at these points. 



3. Bilayer graphene 

We now turn to bilayer graphene and use again the nearest-neighbour, tight-binding 
Hamiltonian in the continuum approximation with k close to the K point. If we 
include a potential difference between the two layers, the Hamiltonian is given by 



n 





Ui 







{ 













VpTT 
U2 



\ 



(3.1) 



Here Ui and U2 are the potentials on layers 1 and 2, respectively, 2 A = Ui — U2 
is the potential difference, and t±_ describes the coupling between the layers. The 



energy spectrum for free electrons is given by ( McCann 2006 Barbier ei a/. 20096) 

+2 

(4A2fc2 



"0 ± 



t 



a/2 



1/2 



t 



i - {4A^k^ + k^t\ 



4 ' 



1/2 



(3.2) 



A. Contrary to Sec. |3]we use units in inverse 
= Uj/hvp, and k = [X^ + kyY^^. This spectrum 
exhibits an energy gap that for 2A <C t± equals the difference 2A between the 



with ui = uo + A and 1*2 — 
distance, namely, e = E/hvp. 



uo 

U4 



conduction and valence band at the K point ( McCann 2006 ) . 

Solutions for this Hamiltonian are four-vectors ij: and for ID potentials we can 
write ip{x,y) — tpi^) exp{ikyy). If the potentials Ui and U2 do not vary in space, 
these solutions are of the form 



*±(a;) 



f± 
h± 
\9±h±J 



(3.3) 
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with /± = [-iky ± \]/[e' - 5], h± = [(e' - 5f - ^ >^^]/[t±{£' ~ S)], and g± 
[iky ± A]/[e' + 6]; the wave vector A is given by 



A+ = 



e'2 + S^- ky^ ± + t\[e'^ - 5^) 



1/2 



(3.4) 



We will write A+ = a and A_ = (3. 



(a) Tuning of the hand offsets 



It was shown before that using a ID biasing, indicated in Figs. 11 |a,b,c) by 
2 A, one can create three types of heterostructures in graphene ( Dragoman et al\ 



2010). A fourth type, where the energy gap is spatially kept constant but the bias 
periodically changes sign along the interfaces, can be introduced (see Fig. [TTJd)). 
We characterize these heterostructures as follows: 

1) Type I: The gate bias applied in the barrier regions is larger than in the well 
regions. 

2) Type II: The gaps, not necessarily equal, are shifted in energy but they have 
an overlap as shown. 

3) Type III: The gaps, not necessarily equal, are shifted in energy and have no 
overlap. 

4) Type IV: The bias changes sign between successive barriers and wells but its 
magnitude remains constant. 



Type IV structures have been shown to localize the wave function at the inter- 



faces ( Martin et al. 2008 Martinez et al. 2009 1 . To understand the influence of 



such interfaces in this section we will separately investigate structures with such a 
single interface embedded by an anti-symmetric potential. 




Figure 11. Four different types of band alignments in bilayer graphene. E^.b, Ec,w, E^^^, 
and denote the energies of the conduction (c) and valence (v) bands in the barrier 
(b) and well (w) regions. The corresponding gap is, respectively, 2Ab and 2Au,. 
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To describe the transmission and bound states of some simple structures we 
notice that in the energy region of interest, i.e., for \E\ < t±, the eigenstates which 
are propagating are the ones with A = a. Accordingly, from now on we will assume 
that /? is complex. In this way we can simply use the transfer-matrix approach of 
Sec. [2] in the transmission calculations. This leads to the relation 












= N 


r 
















(3.5) 



Again the transmission is given by T = 

For a single barrier the transmission in bilayer graphene is given by a com- 
plicated expression. Therefore, we will first look at a few limiting cases. First we 
assume a zero bias A = that corresponds to a particular case of type III het- 
erostructures. In this case we slightly change the definition of the wave vectors: for 
A = we assume a{fi) — [e^ + {—)etj_ — ky]^^^. If we restrict the motion along the 
X axis, by taking ky = 0, and assume a bias A = 0, then the transmission T — |tp 
is given via 



l/t = 



iaoD 



cos{abD) — iQ sm{abD)] 



1 



aoEb 



(3.6) 



This expression depends only on the propagating wave vector a (/3 for E < 0) 
as propagating and localized states are decoupled in this approximation. This also 
means that one does not find any resonances in the transmission for energies in 
the barrier region, i.e., for < e < it. Due to the coupling for nonzero ky with 
the localized states, resonances in the transmission will occur (see Fig. 12). We can 
easily generalize this expression to account for the double barrier case under the 
same assumptions. With an inter-barrier distance one obtains the transmission 
dBarbier et al||2009&|) Td = \td\^ from 



I _ |^|2gi2</)rgi2QoW„ ' 



(3.7) 



with r = |r|e"^' , and t = |t|e"^', corresponding to the single barrier transmission 
and refiection amplitudes. In this case we do have resonances due to the well states; 
they occur for e'^'^'^e'^""^™ = 1. As 0^ is independent of W^^, one obtains more 
resonances by increasing Ww 

For a single i5-function barrier with potential V{x)/hvF = P5{x) under zero 
bias, we find the transmission amplitude 



1/i = cos P -|- i/z sin P - 



[a-pfkl sinP 
Aafie^ cot P + 



(3.8) 



where fi — {e + 1/2) /a and v ~ {e — 1/2)//?. Notice that this formula is periodic in 
the strength of the barrier P as in the single-layer case. 

For the general case we obtained numerical results for the transmission through 
various types of single and double barrier structures; they are shown in Fig. |13| 
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Figure 12. (a) Contour plot of the transmission for the potential of Fig. [ijb) in bilayer 
graphene with Wb = Ww = 40 nm, T4 = —Vw = 100 meV and zero bias. Bound states are 
shown by the red curves, (b) Spectrum for a SL whose unit is the potential structure of 
Fig. [l][b). Blue and red curves show, respectively, the = and kx =-k/L results which 
dehmit the energy bands (green coloured regions). 



The difFcrcnt types of structures clearly lead to different behaviour of the tunnelling 
resonances. 

An interesting structure to study is the fourth type of SLs shown in FigJTT|d) 



To investigate the influence of the localized states (Martin et al. 2008 Martinez 



et al. 2009 1 on the transport properties we embed the anti-symmetric potential 



profile in a structure with unbiased layers. 

Conductance At zero temperature G can be calculated from the transmission 
using Eq. (13) with Go = (4e^Ly/27r/i) {Ep + t^EpY^'^ /hvp for bilayer graphene 
and Ly the width of the sample. The angle of incidence (j) is given by tan0 = ky/a 
with a the wave vector outside the barrier. Figure [T4| shows G for the four SL types. 
Notice the clear differences in 1) the onset of the conductance and 2) the number 
and amplitude of the oscillations. 

Bound states. To describe bound states we assume that there are no propagat- 
ing states, i.e., a and /3 are imaginary or complex (the latter case can be solved 
separately), and only the eigenstates with exponentially decaying behaviour are 
nonzero leading to the relation 







--N 







(3.9) 



From this relation we can find the dispersion relation for the bound states. 

To study the localized states for the anti-symmetric potential profile ( [Martin 



et al. 2008 Martinez et al. 20091 we will use a sharp kink profile (step function). 



The spectrum found by the method above is shown in Fig. [I5[a). We see that there 
are two bound states, both with negative group velocity Vy cx de/dky, as found 



previously by (Martin et al. 2008). No bound state near zero energy was found for 
ky oo in contradiction with ( Martinez et al. 2009 ) . For zero energy we find the 
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Figure 13. (Color online) Contour plot of the transmission through a single barrier in (a) 
and (b), for width Wb = 50 nm, and through double barriers in (c), (d), (e), and (f) of 
equal widths Wb = 20 nm that are separated by Ww ~ 20 nm. Other parameters are as 
follows: (a) = 100 meV, Vt = meV. (b) At = 20 meV, 14 = 50 meV. (c) Type I: 
Vi, = K, = meV, A,„ = 20 meV, and At = 100 meV. (d) Type 11: Vb = -Vw = 20 
meV, A„ = A = 50 meV, (e) Type III: Vb = -K, = 50 meV, A„ = Ab = 20 meV. (f) 
Type IV: 14 = K, = meV, Ab = -A^ = 100 meV. 



solution 



« ±yAiI/23/4, A « i^; 



(3.10) 



the approximation on the second line leads to the expression found by (Martin 



et al. 20081. 



(6) Superlattices 
The heterostructures above (see Fig, 



11), can be used to create four different 



types of SLs (Dragoman et al. 2010). We will especially focus on type IV and type 
III SLs in certain limiting cases. 
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Figure 14. (Color online) Two-terminal conductance of four equally spaced barriers vs 
energy for Wb ~ Ww = 10 nm and different SL types I-IV. The solid red curve (type I) is 
for At = 50 meV, A^, — 20 meV, and Vm = Vb = 0. The blue dashed curve (type 11) is 
for Ai, = All, = 50 meV and Vt — —Vw = 20 meV. The green dotted curve (type 111) is 
for At = Ai„ = 20 meV and Vb = ~Vw ~ 50 meV. The black dash-dotted curve (type IV) 
is for Afc — —Am — 50 meV and Vw ~ Vb = 0. 




Figure 15. (a) Bound states of the anti-symmetric potential profile (type IV) with bias 
A„ = — Ai, — 200 meV. (b) Contour plot of the transmission through a 20 nm wide barrier 
consisting of two regions with opposite biases A — ±100 meV. 



For a type I SL we see in Fig.|l6[a) that the conduction and valence band of the 
bilayer structure are qualitatively similar to those in the presence of a uniform bias. 



Type II structures maintain this gap, see Fig. 16 ^b), as there is a range in energy 
for which there is a gap in the SL potential in the barrier and well regions. In type 
III structures we have two interesting features, which can close the gap. First we 



see from Fig. 12 b) that for zero bias, similar to single-layer graphene, extra Dirac 
points appear for = 0, likewise for Fig. Qd). In the case Wb = Ww = L/2 = W, 
kx = and E = the values for the ky where extra Dirac points occur are given 
by the following transcendental equation 

[cos(aVr) cos(/3M^) - 1] + " i ~ sin(aiy) sm{l3W) = 0. (3.11) 

zap 

Comparing the figures ^h) andgd) we remark that, different from the single- 
layer case, for bilayer graphene the bands in the barrier region are not only flat 
in the x direction for large ky values but also for small ky. The latter corresponds 
to the zero transmission value inside the barrier region for tunneling through a 
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Figure 16. (Color online) Lowest conduction and highest valence band ol the spectrum 
for a square SL with period L = 20 nm and Wb = Ww ~ 10 nm. (a) Type I: At = 100 
meV and A,„ = 0. (b) Type IT. As in (a) for A;, = A„ = 50 meV, and Vt = -14, = 25 
meV. (c) Type lU: Vb = -V^ = 25 meV, and At = = 0. (d) Type IB: 14 = -V^ = 50 
meV and A;, = A^, = 0. (e) Type IV: Plot of the spectrum for a square SL with average 
potential Vb = = and A;, — — A„ = 100 meV. The contours are for the conduction 
band and show that the dispersion is almost flat in the x direction. 



single unbiased barrier in bilayer graphene. Secondly, if there are no extra Dirac 
points (small parameter uL) for certain SL parameters, the gap closes at two points 
at the Fermi-level for ky = 0. The latter we will investigate a bit more in the 
extended Kronig- Penney model. Periodically changing the sign of the bias (type IV) 
introduces a splitting of the charge neutrality point along the ky axis; this agrees 



with what was found by ( Martin et al. 2008 1 . We illustrate that in Fig. [I3[e) for a 



SL with Ah = — A^, = 100 meV. We also see that the two valleys in the spectrum are 
rather flat in the x direction. Upon increasing the parameter Ai, the two touching 
points shift to larger ±.ky and the valleys become flatter in the x direction. For all 
four types of SLs the spectrum is anisotropic and results in very different velocities 
along the x and y directions. 

Extended Kronig-Penney model. To understand which SL parameters lead to 
the creation of a gap we look at the Kronig-Penney limit of type III SLs for zero 
bias (Barbier et al. 2010, unpublished work). Also we choose the extended Kronig- 
Penney model to ensure spectra symmetric with respect to the zero-energy value, 
such that the zero-energy solutions can be traced down more easily. If the latter zero 
modes exist, there is no gap. To simplify the calculations we restrict the spectrum 
to that for ky = 0. This assumption is certainly not valid if the parameter uL 
is large because in that case we expect extra Dirac points (not in the KP limit) 
to appear that will close the gap. The spectrum for ky — is determined by the 
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transcendental equations 

cosfc^L = cos aL cos^ P + sin^ P, 
cos kxL = cos f3Lcos^ P + Dp sin^ P, 



(3.12a) 
(3.12b) 



with Dx = [(A^ + e^) cos AL - + e^] /4A^e^, and X = a, (3. To see whether there 
is a gap in the spectrum we look for a solution with e = in the dispersion relations. 
This gives two values for kx where zero energy solutions occur 



kx.Q = ± arccos[l - {L'^/8) sin^ P]/L, 



(3.13) 



and the crossing points are at (e, kx, ky) = (0, -iikx^, 0). If the kxfi value is not 
real; then t here is no solution at zero energy and a gap arises in the spectrum. From 
Eq. (3.12al we see that for sin^ P > a band gap arises. 



Conductivity. In bilayer graphene the diffusive dc conductivity, given by Eq. ( 2.27 1 , 
takes the form 



ap^(£F)/fTo = (fc|/4^4) [l ± S/2{klS + 1/4)1/2 



(3.14) 



with kp = [el + A^T {spS ~ A2)i/2]i/2^ ,5 = 1 + AA^, and ctq = e^Tpt^/h^ 




Figure 17. (Color online) Conductivities, a^x in (a) and ayy in (b), vs Fermi energy for 
the four types of SLs with L = 20 nm and Wt = Ww = 10 nm, at temperature T = 45K; 
(70 = e^Tpt^/h^. Type I: Ab = 50 meV, A„ = 25 meV and Vt = V^u = 0. Type U: 
Ab = A» = 25 meV and Vb = -V^ = 50 meV. Type HI: At = A„ = 50 meV and 
Vb = -Vu, = 25 meV. Type IV: Ab = -A„ = 100 meV and Vb = K„ = 0. 



In Figs. 17 |a), (b) the conductivities a^x in (a) and ayy in 



graphene are shown for the various types of SLs defined in SecKMiy 



(b) for bilayer 
Notice that for 



type IV SL the conductivities ax 
in the spectrum. 



and ayy differ substantially due to the anisotropy 



4. Conclusions 

We reviewed the electronic band structure of single-layer and bilayer graphene in 
the presence of ID periodic potentials. In addition, we investigated the conditions 
that lead to carrier collimation in single-layer graphene and determined when extra 
Dirac points appear in the spectrum and what their influence is on the conductivity. 
Furthermore, we investigated the tunnelling through, and bound states created by. 
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simple barrier structures. In single-layer graphene we found that the SL spectrum 
can be linked to the bound states of a combined barrier and a well. 

In bilayer graphene we considered transport through different types of het- 
erostructures, where we distinguished between four types of band alignments. We 
also connected the bound states in an anti-symmetric potential (type IV) with the 
transmission through such a potential barrier. Furthermore, we investigated the 
same four types of band alignments in SLs. The differences between the four types 
of SLs are reflected not only in the spectrum but also in the conductivities parallel 
and perpendicular to the SL direction. For type III SLs, which have a zero bias, 
we found a feature in the spectrum similar to the extra Dirac points found for 
single-layer graphene. Also, for not too large strengths of the SL barriers we found 
that the valence and condunction bands touch at points in k space with ky — and 
nonzero ky. Type IV SLs tend to split the K (K') valley into two valleys. 

In the Kronig-Penney limit, where we take the barriers to be S functions V{x) /hvp = 
PS{x), we saw that the SL spectra, the transmission, conductance, etc., are periodic 
in the strength of the barriers. As s well known, this is not the case for standard 
electrons. An important qualitatively new feature is encountered in the extended 
Kronig-Penney limit for P = (n + l/2)7r, see Sec. [2j|d]): the Dirac point becomes a 
Dirac line. 

We expect that these relatively recent findings, that we reviewed in this work, 
will be tested experimentally in the near future. 

This work was supported by IMEC, the Flemish Science Foundation (FWO-Vl), the Bel- 
gian Science Policy (lAP), and the Canadian NSERC Grant No. OGP0121756. 
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